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ABSTRACT 



The purpose of this report is to present a breakdown surface model for 
evaluating thermal backscatter ing flux from the supersonic exhaust plume of a 



space-based HF laser. 


The plume is of ring symmetry. It consists of .a 


gaseous mixture of H, 


HF, H2 , DF and He. Fluxes of these species are 


considered separately. 


The model is carefully analyzed and is shown to 


overestimate the flux. 


Actual flux levels of the heavy corrosive molecules 



(HF, DF) have been found to be exceedingly low. It is concluded that the 
contribution of thermal backscattering to contaminating flux of HF and DF can 
be neglected. This work is an extension and modification of the recent thesis 



work done by S. E. McCarty at the Naval Postgraduate School 
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1. INTRODUCTION 



This report is a presentation of one part of a study on the contaminating 
backflow from the exhaust plume of a large space-based HF laser (Figure 1). 
The exhaust plume is an underexpanded supersonic ring-jet, designed to stay 
clear of the spacecraft by maintaining a Prandt 1-Meyer turning angle at the 
nozzle lips of well below 90°. However, it is wpll known from experience with 
rocket, plumes in spacell ,2] that cavity regions (where continuum gasdynamic 
theory predicts vacuum) are filled with a free-molecular flow. This back flow 
is largely due to viscous effects, which give rise to a M spill-over M of the 
boundary layer around the nozzle lip[3]. Assuming the boundary layer can be 

eliminated (e.g., an expanding step design of the nozzle lip), there are 
two more mechanisms which lead to backflow: thermal backscattering and 

scattering by ambient molecules traveling at orbital speeds* Since these 
effects are a small perturbation to the exhaust flow field, they can be 
considered independently (the total backflow will be a superposition of 
contributions due to these two effects). As a first phase of our broader 
study, we consider solely the contribution of thermal backscattering to back- 
flow from a ring-plume of an HF laser, via a simple model of molecular 
effusion from a breakdown surface, fashioned after ideas suggested by 
NollerlM. Our results indicate that the backflow of the heavier contaminants 
(HF, DF) due to thermal backscattering is negligible. 

Naturally, our study pertains to presumably typical operating conditions 
of the HF laser. These operating conditions were largely determined from a 
report on some HF laser tests conducted at TRW in 1971 [5] (in particular, 
Table 5, Test III, of this report). The typical parameters at the nozzle^exit 
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are : 

Composition (mole fractions): 

Specific Heats Ratio: 

Mach Number 

Average Molecular Weight 

Stagnation Temperature 

Stagnation Density 

Molecular Diameter, assuming it is 
uniform for all species (hard- 
sphere collisions ) 

The exit Mach number can be chosen higher than Mj * 4, but not 
considerably lower than this value, since M^ * 4 results in a modest clearance 
angle of 41° between the limiting (vacuum) characteristic of the lip-centered 
rarefaction fan and the spacecraft. We assume isentropic flow throughout the 
diffuser^, so that upon specifying the composition and flow variables at the 
diffuser inlet, along with M^ at the diffuser exit, the exit flow is fully 
determined. One exception to this definition, however, is the stagnation 
temperature, which was estimated as T 0 - 1400 [K] at the diffuser inlet 
We set T 0 = 2300 [K], which corresponds to complete hydrogen recombination, 

even though the flow in the diffuser is of a nearly frozen composition due to 
the low rate of hydrogen recombination^]. The reason for this choice is that 
given the uncertainty in determining T 0 , which results from an uncertainty in 
the degree of hydrogen recombination, it is the most pessimistic choice, 
resulting in higher thermally backscattered flux. 

*> M 

The model that we propose for evaluating the backscattered flux arriving 



[H] = .091, [HF] - .091, 

[H 2 ] = .104, [DF] = .135, 

[He] « .579 

Y = 1.54 (assuming ideal gas) 
Hi = 4.0 

W A = 7.27 [kg/kg mole] (1.1) 
T 0 - 2300 [K] 

P 0 = 0.0075 [kg/m 3 ] 

D = 2.5 x 10 -10 [m ] 
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at the spacecraft (Figures 1, 2) is based on the effusive breakdown surface 
concept suggested by Nollerl^l. The gradual transition from continuum to 
collisionless flow, which invariably takes place at the outer fringes of 
exhaust plumes having a near-vacuum background environment, is replaced by an 
abrupt change* We assume that the flow regime in each stream tube changes 
from continuum (with local thermodynnamic equilibrium) to collisionless, upon 
crossing some breakdown surface* 

An important simplification is introduced in the case of a large-radius 
spacecraft (about 2*5[m]), by observing that the temperature along the 
breakdown surface decreases so- sharply with the distance from the . nozzle lip, 
that the segment contributing significantly to thermal backscattering is only 
about 0.01 to 0.1 [m] long. Consequently, the lip-centered rarefaction 
ring-fan may well be approximated by the standard (planar) Prandt 1-Meyer 
flow field. 

The structure of this report is as follows. The breakdown surface and the 
molecular effusion flux from it are obtained in closed-form expressions in 
Section 2. Section 3 is devoted to the spatial integration scheme, which is 
the evaluation of the flux arriving at a certain point on the spacecraft. 
Results of flux distribution along the spacecraft for the presumed laser 
operating range are presented and discussed in Section 4, followed by a 
critical examination of the breakdown surface, model. Conclusions are given 
in Section 5, and Section 6 is a list of references. The code RINGBD , which 
computes the flux by numerical integration over the breakdown (effusing) 
surface, is given in Appendix A. 
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2 . BREAKDOWN SURFACE AND EFFUSION FLUX 



Our model for the thermally backscattered flux arriving at the surface of 
the spacecraft is essentially a modification of Noller f s concept of a 
breakdown effusive surface^. We substitute his definition of a breakdown 
surface by a similar one introduced by Bird^> Section 8.3] # We obtain the 

one-sided effusion flux from the breakdown surface by integrating over 
velocity space as suggested by Noller^, except for the fact that we compute 
flux rather than density and we also consider the flux of species having 
molecular weight different from the average. In the following, each one of 
these steps is described in some detail, beginning with the breakdown surface. 

As mentioned in the introduction, the lip-centered rarefaction fan is 
approximated by a planar Prand tl-Meyer flow field (Figure 2). The standard 
expressions for this flow field have Mach number (M) as the independent 
parameter, thus M varies between M = Mj at the exit and M + ® at the limiting 
(vacuum) characteristic. (Index 1 always refers to exit conditions, i.e., to 
parameters evaluated at M = M^). 

*(m) = - c(m) + C 1 + 1^+ “ 

C(M) = T 1/2 ARCTAN[r -1/2 (M 2 -l) 1/2 ] r= (^Tf) 

u(M) = ARCS IN( M -1 ) (2.1) 

0(M) = \p(M) - ii(tl) ' 

where is the angle of characteristic lines, and 0 is the angle of the 
velocity vector (or streamline). 
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Adopting Bird's definition of a breakdown parameter, which was first 
introduced in conjunction with a spherical source flow^ > Section 8.3] anc j 
later was shown to be meaningful also in a Prand tl-Meyer flow^] , we define 
the breakdown surface as having a constant value of *B, where B is given by: 



= U | dp I 
v p I dS I 



( 2 . 2 ) 



Here p, U, v, S are local flow density, speed, collision frequency, and 
coordinate along streamlines (thus restricting this definition of B to 
stationary flows). From the geometrical relationships in a Prand tl-Meyer fan 
(Figure 2) and from (2.1) we get: 

If = - (1/R) (If) sin » = - t M ' 1(fl2 - 1)1/2 1 ( p /R) 

1 (2.3) 

P(M) = Po (1 M2 ) 

Using the expression for collis ion frequency ^ 6 ] ; 



V 0 = 4(«/y) 1/2 (N 0 D 2 C q ) (2.4) 

where N 0 ,C Q ,D are stagnation number density, stagnation sound speed, molecular 
diameter, and using U = MC in conjunction with (2.2) and (2.3), we get: 

1/2 _L 

R b (M) = (BN 0 D 2 )-1 (M 2 -l ) 1/2 (l+^y- M 2 ) Y-l (2.5) 

This expression is almost identical to that of Bird 1 , the main difference 
being in assuming a constant collision diameter (hard spheres), which we 
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believe to be commensurate with the overall crudeness of the model. The 



breakdown surface as defined by equations 2,5 and 2.1, starts at point 
[RgCM^), \Jj(M ^ ) ] on the exit characteristic M = , which we refer to as the 

intial point (see Figure 4). However, a breakdown- in continuum flow also 
takes place on the segment of the exit characteristic between the corner and 
the intial point, since the value of the breakdown parameter there (Equation 
2.2) is clearly larger than the value of B ufeed in defining the breakdown 
surface (Equation 2.5). Hence, the breakdown surface defined by (2.5) should 
be supplemented by that segment. We refer to the combined surface as the 

augmented breakdown surface. The segment on the exit characteristic is 
referred to as the supplementary breakdown surface. 

The one-sided directed effusion flux is defined as the number flux of 
molecules per unit area of an area element normal to the flux direction, per 
unit solid angle about the flux direction. It is obtained as a function of 
local Mach number and the angle < between the flux direction and the local 
velocity vector, by repeating Noller's velocity integration scheme^» E Q* (6)], 
with an added factor of molecular speed in order to obtain flux (rather than 
density as in Noller's work). The resulting expression for species i is 

readily obtained by using standard definite integrals: 

F t (M) = h* N 0 C q (W A /W i ) 3 ^ 2 (! + lZk M 2 ) 2(y_1) 

[(2yTt 3 )" 1/2 (l + ( 1 / 2) ytl 2 cos 2 < ) EXP (- yM 2 /2) + 

( 2 . 6 ) 

( 2tt )"* 1 M (3/2 + (l/2)yM 2 cos 2 <) cos< EXP (- (l/2)yM 2 sin 2 >c) 

1/2 ~ 

ERFC (- (y/2) M cosk)] 

00 

ERFC(v) = 2 1/2 | £xp (- x 2)^ x (complementary error function) 

v 
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(2.6 Continued) 



M = (Wi/W A ) 1/2 M 

- Hole fraction of species i. 

- Molecular weight of species i. 

The dependence of Fj/M) on the flux angle < is so sensitive that for some 
Mach number around M « 10* the backflow (in the typical operating range) is 
virtually negligible* In the following section we describe how the flux F^(M) 
is integrated over the augmented breakdown surface, yielding the backscattered 
flux arriving at the surface of the spacecraft. 
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3 . FLUX INTEGRATION 



The effusion flux F^(M) given by (2.6) above, is defined in such a way 
that the number of molecules effusing from an area element AAg of the break- 
down surface and arriving at an area element AAg on the spacecraft (per 
second), is given by: 

F^(M) ( AAg COS ag) (AAg COSag) L ^ [molecules per second] (3.1) 

D o 

where a g , a g are the angles between the line-of-sight Lgg (Figure 3) and the 
normals to the breakdown surface and the spacecraft surface respectively. Lgg 
is the distance between the elements AAg and AAg. Dividing equation (3.1) by 
AAg and integrating over the breakdown surface, the flux per unit area of 
the spacecraft is given by: 

Qi = |F^(M)COS a g COS a g L ^ dAg [Molecules per second per m^] (3.2) 

3 S 

The integration scheme for over the breakdown surface is expressed in 
terms of the set of polar coordinates R, \J;, $ (Figure 3). For a point 

(R, i|i , <j>) on the breakdown' surface, using Cartesian coordinates (X, Y, Z) and 
the angle u between X-axis and the line-of-sight Lgg we obtain the following 
geometrical relationships: 

X = R COS 4; ; Y = (A 0 + R SIN<|>) COS^j ; Z = (A Q + R SINi|>)SIN<f>- 

A o 

C0S<j, MAX (H) = ^Ao+R(M)SIN 4») (3.3) 

U = U(COS0, SIN0 COS<|>, SIN0 SINt})) 

Lbs/I^BsI = (COSu, SINoj COSB, SINcj SIN 0) 

TANB 
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The cosines COSk, COS a g, COS a g in ( 3.2 ) are expressed as scalar produces 
of (Lnc/lLssN and unit vectors along the local velocity vector U, along 
the local normal to the breakdown surface and along the local normal to the 
spacecraft surface, correspondingly. 

The integration is performed .numerically in two phases, the first being 
the integration along the supplimentary breakdown surface (Figure 4), For 
this first phase, the straight line segment which constitutes the 
supplementary breakdown surface is divided into several intervals of length AR 
(typically 10 intervals). Each interval generates a half-strip by rotating it 
from <j> = 0 through § = This strip is in turn subdivided into 
several sub-intervals of A<() each (typically 10 intervals). The total flux 
arriving at X s is obtained by summing contributions from each sub-interval 
(two-dimensional integration). When the integration along the supplementary 
breakdown surface is concluded, it is continued into the breakdown surface, 
where AR intervals are replaced by breakdown surface intervals that correspond 
to a fixed Mach number increment AM (typically AM=0.1). The integration 
proceeds along the breakdown surface (2.5) until the contribution of the last 
AM strip is negligibly small. The computation time is modest (about 1 second 
CPU per X s point, on IBM 3033 mainframe computer). The computations were 
carried out by a code RINGBD written specifically for this purpose. Further 
details of the scheme and programming can be obtained by reading this code 
which is given in Appendix A. 
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4. RESULTS AND DISCUSSION 



4.1 Presentation of Results 

The molecular flux backscattered to the spacecraft from the surface 
of continuum flow breakdown in the lip-centered rarefaction fan, has been 
computed for all five species H, HF, H 2 > DF, He. The results are depicted in 
Figures 5 to 9 respectively. For each species two more cases were computed in 
addition to the nominal case (1.1), where the stagnation density p 0 was 
replaced by p o /10 and by p o *10 (see Figures 5 to 9). This has been done in 
order to demonstrate the effect of variations in exit flow conditions on the 
flux. The particular choice of p 0 was motivated by the fact that the effects 
of changing p 0 are not obvious. The effects of changing the exit Mach number 
M or the stagnation temperature T 0 are rather obvious (a higher flux would 
result from either a decrease in M or an increase in T 0 ) • It turns out that 
for points lying not too near the nozzle lip (X s > 0.1 m) , the lower density 
flow generates a higher backscattered flux! 

In addition to varying p 0 , we also varied the breakdown parameter B, 
obtaining a surprising result. The computation was performed for a particular 
species (HF), and the results obtained upon replacing B=0.05 (nominal value) 
by B/2 and by B*2 are brought in Figure 10. 

It turns out that the B/2 case has the higher flux. This is somewhat 
surprising, since a lower value of B in a centered rarefaction fan (equation 
2.5) means that the breakdown of continuum flow takes place in a region 
further out from the corner. In a source flow (e.g., a spherical source), 
that implies lower density and temperature, which would give rise to 
lower thermally backscattered flux. 
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An explanation to these seemingly counterintuitive results, along 
with some deeper insight into the breakdown surface model as it is applied to 
a centered rarefaction flow, can be obtained by taking a close look at the 



flow field and the breakdown surface in the vicinity of the corner. We take 
up this matter in the following sections. 



the nominal case) of the five species with each other (Figure 11). This 

figure underlines the fact that the flux of light species (H, H 2 > He) is many 
orders of magnitude (typically 10 15 ) times that of heavy species (HF, DF) . 

Indeed, these results demonstrate a well known effect: When an expanding 

gaseous mixture of light and heavy molecules experiences a breakdown of 
continuum flow, a separation of species takes place (see e.g., the work of 
Cattolica et. al. [8]). 

4.2 The Breakdown Surface and Streamlines 

Consider the parametric description Rg(M) for the breakdown surface 
(Equation 2.5). Normalizing R relative to the exit mean free path Xj , we get: 



We conclude the presentation of results, by comparing the flux (in 



1 



R B CM> - Rb(M,)[(m2 - 1)/(M, - 1)] 1/2 [(1 + J=i m2)/{1 + M 2 )]-'- 1 

R B Oii)/x, - [(Y./2) 1/2 /{r«)B] (m 2 - i) 1 ' 2 



(4.1) 



1 




The normalized surface Rg(M)/X^ is thus independent of stagnation density, 
depending only on y, M 1 , B. 
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Let us now derive a parametric equation Rg(M) for a streamline that 
enters the fan at point R^M^) on the exit characteristic. The following 
geometrical relationship is readily obtained by considering two characteristic 
lines and and a streamline inclined at the Mach angle p to them: 



dRgC^) 

dtp 



= - Rg(ij>) (tanu)-l 



Using the standard Prandtl-Meyer functions (2.1), we 
following differential equation for Rg(M): 



(4.2) 
get the 



1 

r^TmT 



M(1 



2.1 

2 



M 2 ) 



_1 



(4.3) 



This equation is readily integrated, giving: 

. R S (11) = Rg ( ) [(l + M 2 )/(l + 



fVn 



(4.4) 



As pointed out by Birdt^l, there is a particular streamline R sa (M) 
which asymptotically approaches the breakdown surface for large M, since the 
ratio RgOO/RgCM) tends to a constant (not zero) when M «. (Strictly 
speaking, this holds only for hard-sphere molecules, i.e., only when w = 0.5 
in 1 ) . The limit is : 



lim Rs(M) Rs(Mi> f M? -1 P /2 

r^TmT = + (4 ' 5) 

1 y-1 

For the limiting streamline Rg a (M) , the ratio R sa (M)/Rg(M) should 
tend to 1. This determines the point Rg a (M ) at which the limiting streamline 
enters the fan, as well as the entire line Rg a (M) : 






R B (M 1 ) 



2 

-Ml 



Ilf - i 



1/2 



(4.6) 
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Clearly, R sa (M) is larger than Rg(M) for any H > Mp so that no 
streamline beyond R sa (M) can cross the breakdown surface. This pattern is 
shown in Figure 4, where Rg a (M) is denoted "streamline 2", and "streamline 1" 
is the streamline R^M^ = RgCM^. 

All this leads to the following observation regarding the continuum 
breakdown of the flow in a centered rarefaction f an ^ • Referring to Figure 
4, the fluid entering the fan through the supplementary breakdown surface 
(i.e., through the exit characteristic between the corner and the initial 
point of streamline 1), experiences breakdown immediately upon crossing this 
surface. Every streamline between streamline 1 and streamline 2 crosses the 
breakdown surface at some Mach number M > M , and at that, point the 
continuum flow breaks down. All fluid entering the fan beyond stream- line 2 
will never pass through the breakdown surface, and hence will maintain a 
continuum flow regime all the way to infinity. Of course, that is only true 
for planar centered rarefaction fans. When the exhaust flow emerges from a 
nozzle of finite width, and especially when the exhaust jet has a ring 
symmetry (as in our case), the breakdown surface gradually curves in a 
balloonlike shape towards the opposite nozzle lip, forming the familiar plume 
pattern (Figure 1). 

4 . 3 Ana lysis and Discussion of Results 

The foregoing analysis is now used to explain the variation in back- 
scattered flux due to a change in exhaust flow conditions at the nozzle exit. 
Specifically, we consider a tenfold decrease in stagnation density* (i.e., the 
case p o /10), and hence a tenfold increase in the exit mean free path X^. 

The effusion flux from the breakdown surface is proportional to the 
local density, so one would expect to observe a decrease in flux, rather than 



13 



an increase (see Figures 5 to 9 , for X s >0.1 m) . Other factors causing 
increased flux, must then be larger than 10 so that they more than offset the 
1/10 factor in density. It turns out that these effects are mainly 
geometrical, in that a tenfold increase in X^ causes the domain of integration 
on the breakdown surface to increase more than tenfold. In the (X,Y) plane 
there is a tenfold "blowup** of the breakdown surface, due to the self-similar 
structure of the Prand tl-Meyer flow field. As a result of this ‘’blowup" in 
(X ,Y) , the angular integration range $> max also increases, albeit not linearly 
(Equation 3.3). Another geometrical effect is an increase in the flux 
incidence cosine factor cosa s (see Equation 3.2), which for points X g 
sufficiently far from the nozzle lip, increases roughly tenfold (while the 
other cosine factor cos a g is almost constant). All this provides a* 
qualitative explanation for the observed increase in flux at far points (X s > 
0.1 m) . 



As for the near range (X s < 0.1 m) , another effect becomes 
increasingly significant as X s approaches the nozzle lip. The turning angle 
k, by which backseat tered molecules have to be deflected relative to the flow 
velocity vector in order to reach point X s on the spacecraft (Figure 2), 
increases with the size of the breakdown surface (fixed M* and X s ). Since the 
local effusion flux (Equation 2.6) decreases rather sharply as k is increased, 
the net result is a tendency to get a reduced backscattered flux at near 
points such as Xg = 0.01 m (Figures 3 to 9) . 

We now turn to the effect of changing the value of the breakdown 
parameter B. From equation 4.1 it is clear that multiplying B by some factor 
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will have the same “blowup" effect as dividing M by the same factor, A 

tenfold decrease in B is thus geometrically equivalent to a tenfold decrease 

in p 0 * However, since the local effusion flux at the breakdown surface is 
proportional to p Q while it is independent of B, the B/10 case will have ten 
times as much backscattered flux as the p o /10 case. In order to illustrate 
the sensitivity of the flux estimates to an uncertainty in the appropriate 
value of B, we computed the cases B/ 2 and B*2 for one species (HF), and they 
are presented in Figure 10. The variation in flux relative to the nominal 
case (B = 0.05), is by a factor no larger than about 5. Results for other 

species were found to exhibit comparable variations. 

Does this observation about the dependence of the breakdown surface 

on B agree with the breakdown surface appropriate to the far field of the 

exhaust plume? In stationary source flow into vacuum, and when M » 1, the 

5- 1 

breakdown parameter varies with radius as B ~ R (5=1 for cylindrical 
source, 5=2 for spherical source). In a ringjet, the stream tubes of the 
exhaust plume generally diverge at a rate higher than that of stream tubes in 
a cylindrical source flow, so the effective value of 5 in a ringjet is 5 > 1. 
Hence, in this case the far field breakdown surface moves downstream along 
each stream tube as the value of B increases. This is indeed geometrically 
compatible with the fact that near the corner of the lip-centered rarefaction 
fan B ~ R , as shown schematically in Figure 12. The dependence of the break- 
down surface on B near the corner and in the far field, thus assures that 
complete breakdown surfaces corresponding to different values of B, do not 
intersect (Figure 12). 

In the foregoing discussion it was pointed out that variations in 
flux caused by changes in parameters such as p Q and B, were directly related 
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to the self-similar structure of the Prandtl-Meyer flow field 



It has been 



further shown that these variations are well-understood within the framework 
of the breakdown surface model and that they are not excessively large. Are 
we to conclude that the thermally backscattered flux estimates of the present 
model are also physically plausible and reliable? In the following section we 
take up this matter, arriving at some interesting conclusions about this model 
and its range of validity. 
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4.4 Critical Examination of the Model 



Consider the centered rarefaction flow field of a compressible fluid 
negotiating an expansive corner at supersonic speed (Prandt 1-Meyer flow). The 
streamlines of this flow field have an orderly "layered 1 * structure, with each 
streamline curving around the corner, starting at its point of entrance into 
the fan (see Figure 4). 

The present model is based on the stipulation that there is a point 
of continuum flow breakdown on each streamline, provided this streamline is 
not beyond a certain limiting streamline. Consider a sample molecule effusing 
from this breakdown point toward the spacecraft. It advances at constant 
speed along a straight line trajectory, traversing all inner streamlines. 
Since the flow velocity vector points away from the spacecraft, and since the 
flow is highly supersonic so that the velocity of most individual molecules 
does not differ much from the flow velocity (i.e., it is a "cold" flow), any 
collision of the sample molecule with a mainflow molecule will most probably 
divert the sample molecule away from the spacecraft.. What is the probability 
that a sample molecule would traverse this cross flow collisionlessly ? This 
probability is simply exp(-n), where n is the expected number of collisions 
along the straight-line trajectory from the point of breakdown to the 
spacecraft. In the typical operating conditions assumed here, we estimated n 
to be roughly about 10. Since this no-collision probability factor is ignored 
in the formulation of the present model, the backscattered flux may be 
exaggerated by a factor of exp (10) or about 10 4 . We conclude that in all 
likelihood, the prediction of the breakdown surface model for thermally 
backscattered flux from a centered rarefaction flow, is substantially 
overestimated . 
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Can anything be done to improve the present model? One may be 
inclined to suggest at this point that the obvious remedy is to incorporate 
the no-collision probability factor into the model. Rather, we prefer to 
retain the breakdown surface model in its present form as a simple means of 
obtaining an overestimate to the thermally backscattered flux from a centered 
rarefaction flow. An improved model can be constructed by considering thermal 
backscattering from the entire flow field (tempered by the probability of 
no-collision), without resorting to the physically untenable notion of an 
abrupt transition from continuum flow to free molecular flow. 
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5* CONCLUSIONS 



Some surprising similarity laws of the breakdown surface model were 
observed. It has been shown that they were a direct result of the self 
similar structure of the Prandtl-Meyer flow field to which the model was 
applied. Specifically, it was found (and shown plausible) that reduced values 
of either the exhaust stagnation density p 0 , or the breakdown parameter B, 
caused higher backseat tered flux. 

The breakdown surface model for thermally backscattered flux from a 
centered rarefaction fan, has been shown to overestimate the flux arriving at 
the spacecraft. It is suggested that an improved model be constructed by 
considering thermal backscattering from the entire flow field, along with the 
probability factor for a side-scattered molecule traversing the main flow 
collisionlessly . 

The molecular flux of corrosive species (HF, DF) arriving at the 
spacecraft (Figures 6 and 8) is no larger than about 10 7 (sec -1 m“ 2 ), which is 
negligible since it corresponds to about 10 -5 molecular monolayers per year. 
This conclusion is reliable since even this flux level is an overestimate. 

The maximum thermally backscattered flux of light species (H, H2, He) is 
in the range of 10 20 to 10 22 (sec -1 m~ 2 ) (see Figures 5, 7, 9). Thus, we 
conclude that while thermal backscattering would contribute significantly to 
the flux of light molecules arriving at the spacecraft, it is utterly 
negligible as far as heavy molecules are concerned. 
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APPENDIX A. The Computer Code RINGBD 



We present a printout of the code RINGBD along with the reults of the 
nominal case (printout of actual run). This is preceded by a brief 
description of the subroutines and a summary of major variables with their 
code and report notations. 

A. 1 Description of Subroutines 



MAIN PROGRAM 


Computes flux integration by summation of segment 
contributions (centered). Printing of results. 


INI DAT 


Definition of all data (no input file). 
Preparatory evaluation of parameters. Printing 
of data. 


FLUX 


Evaluates flux emitted from a single point on 
breakdown surface (mean segment values) to a 
point on spacecraft (XS). 


BREAKR 


Computes point on breakdown surface for given 
Mach number. 


BREAKM 


Computes point on breakdown surface for mean 
Mach number of a segment. 
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A. 2 Code Versus Report Notation 



XC(I) - [ABj - 

WC(I) - W A 
WAV - W A 
TO - T 0 
RHOO - p 0 
G ~ Y 
EMI - M 1 
LAMDA1- X - 
AO - Aq 
R - R 
DIST - L bs ~ 

XS - Xg - 
PSI - 
AM U - y 
TETA - 0 
PHI - <{> 

W - 03 
BETA - g 

DM0 - AM 
EM - M 
PBIRD - B 



Mole fraction of species AB. (1=1, 2, 3, 4, 5 cor- 

responds to H, HF, H 2 > DF, He). 

Molecular weight of species i. 

Average molecular weight 
Stagnation temperature 
S tagna t ion dens i ty 
Specific heat ratio 
Exit Mach number 
Exit mean free path 
Spacecraft radius 

Distance from corner (X=0, Y 2 +Z 2 = Aq 2 ). 

Distance between emitting point on breakdown surface 
and receiving point (XS) on spacecraft. 

Point on spacecraft (X=X S , Y=A Q , Z=0). 

Characteris tic angle 
Mach angle 

Velocity vector angle 

Rotation angle for flux integration 

*► 

Angle between x-axis and line-of-sight L^g 

Angle between Y-axis and projection of L Bs on (Y,Z) 

plane . 

Mach number increment for flux integration 
Mach number 
Breakdown parame ter 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 



A. 3 Code Listing (Run of Nominal Case) 



$ JOB 



101 

C 

C 

no 



45 



41 

42 



RINGBD, NOXREF 
IMPLICIT REAL*3(A-H,0-Z,$) 

COMMON /GAMA/G,G1,G2,G3,G4,G5,G6,G7,G8,G9,G10,G11,G12,G13,G1 
1 G16,G17,G18,G19,G20 

COMMON /PAR/CO, EN0, EMI, D,TLIM, ETALIM, CL IM, EL0,Q0,T0, 

1 PBI RD, RBI RD, DM0, DEG, OMEGA, XSV( 51) 

COMMON /NPAR/NETA,NC,NT,NEMO,NPHI,NXS,NRO,NSPEC 

COMMON /GE0M/APF,PAI,PAI2,SW,CW,BETA,SBETA,CBETA,PSI1,SPSI1, 

1 CPSI1,PSIF,SPSIF,CPSIF,AK,SK,CK,A0,RF,XF,YF,ZF, 

2 PHI,SPHI,CPHI,RMIN,RMAX,XS,DIST, 

3 AMU1 , ZETA1 , XN , YN , ZN , PSIM , SPSIM, CPSIM, RO 
COMMON /EPSIL/EP5Q, EPSETA, EPST, EPSC, EPSEM 

COMMON /EXT REM/ TEXT , ETA EXT , CEXT , REXT , PSI EXT , EM EXT , BEXT , QEXT 
COMMON /SPEC/NAV,XC(5),WC(5),WCR(5),XNAME(5),QC(5),FLUXC(5) 
DIMENSION DSUMC5) 

PRINT 101 

FORMAT C 1 1 f / IX, 'RINGBD - FLUX INTEGRATION FROM BREAKDOWN 1 , 

1 IX, ’SURFACE 1 //) 

CALL INIDAT 



PRINT 110, XNAME 

FORMAT (///IX, 1 NX 1 ,’ NEM’,’ XS ’, 

1 5C4X,A6,1X, ’/ LOG’, IX)) 

DO 200 NX= 1 , NXS 
EM=EM1 

CALL BREAKRC EM, RF) 

IF(NRO.GT.O) RF=R0 
XF=RF*CPSI1 
YF=RF*SPSI1+A0 
XS=XSV ( NX ) 

QMAX = 0 . 

DO 45 N=1 , NSPEC 
FLUXCC N ) = 0 . 

CONTINUE 

DO 1 NEM=1 , NEMO 

RN = RF 

XN = XF 

YN=YF 

IF(NEM.GT.NRO) GO TO 41 

RF=RO+DFLOAT(NEM)X(RMIN-RO)/DFLOAT(NRO) 

XF=RF*CPSI1 

YF=RF*SPSI1+A0 

RMEAN=(RN+RF)/2.D0 

EMMEAN=EM1 

PSIM=PSI 1 

SPSIM=SPS I 1 

CPSIM=CPSI 1 

GO TO 42 

CONTINUE 

EM=EM+DM0 

CALL BREAKRC EM, RF) 

EMMEAN=EM-DM0/2 . DO 
CALL BREAKMC EMMEAN, RMEAN) 

CONTINUE 

ALONG=DSQRT( (XF-XN)**2+(YF-YN)x*2) 

SAL FA= (YF-YN )/ ALONG 

CA L FA = ( XF-XN )/ALONG 

PHIMAX = DARCOS ( AO/ CAO + RMEANXSPSIM)) 

DPHI=PHIMAX/NPHI 



PHIMAX 1 



QMAX 



RIN00010 
RI NO 0 020 
4 , G15, RIN00030 
RI NOO 040 
• RIN 00050 
RIN00060 
RI NO 0 07 0 
RIN00080 
RIN 00090 
RIN00100 
RI NOO 1 1 0 
RIN00120 
RIN 00130 
RIN00140 
RIN 00150 
RIN00160 
RIN 00170 
RIN00180 
RIN00190 
RIN00200 
RI NO 021 0 
RIN00220 
’, RI NO 0230 
RIN 00240 
RI NO 025 0 
RIN0026 0 
RI NO 027 0 
RIN0028 0 
RI NO 029 0 
RIN 00300 
RIN0031 0 
RIN00320 
RIN 00330 
RIN 00340 
RIN00350 
RIN0036 0 
RIN0037 0 
RIN00380 
RIN 00390 
RIN 00400 
RIN00410 
RIN00420 
RIN00430 
RIN00440 
RIN00450 
RIN0046 0 
RIN 00470 
RI NO 048 0 
RIN00490 
RIN0050 0 
RIN0051 0 
RIN00520 
RIN00530 
RI NO 054 0 
RIN00550 
RIN0056 0 
RIN00570 
RI NO 0580 
RIN 00590 
RIN00600 
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51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

8 4 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 



DO N=1 , NSPEC 
DSUMC N ) = 0 . 

44 CONTINUE 

DO 2 NP = 1 , NPHI 
PHI=CDFLOATCNP)-0 .5D0)*DPHI 
CALL FLUXC EMMEAN, RMEAN ) 

CROSS 1 =SW*CB ETA 

CR0SS2=CSALFA)*C-CW)+C-CALFA*CPHI)*C-SW*CBETA)+ 

1 C -CAL FAXSPHI ) * ( -SW*SBETA ) 

GOREM=CROSS1*CROSS2*DPHIXCAO+RMEAN*SPSIM)*ALONG/DISTX*2 
DO 24 N=1 , NSPEC 
DSUMC N)=DSUMCN)+QCCN)*GOREM 

24 CONTINUE 
IFCQMAX.GE.QEXT) GO TO 25 
QMAX=QEXT 

25 CONTINUE 

C PRINT 22, NEM, NP, EMMEAN, RMEAN, PH I MAX* DEG, DARCOSCCW)*DEG, 

C 1 BETA*DEG,PHI*DEG,CR0SS1,CR0SS2,AL0NG,DIST,G0REM, 

C 2 SPSIM,AO,(QC(N), FLUXCC N ), N = 1 , NSPEC) 

C22 FORMAT (/ IX, 1 NEM, NP, EMMEAN, RMEAN, PHIMAX= ' , 213 , 3D13 . 4/ 

C 1 IX, »W,BETA,PHI=' ,3D15.5/ 

C 2 IX, *CR0SS1,CR0SS2,AL0NG,DIST,G0REM= , ,5D15.5/ 

C 3 IX, , SPSIM,A0= f ,2015.5/ 

C 4 IX, , QC,FLUXC= , ,5(1X,D10.3,1X,D10.3)) 

2 CONTINUE 

DO 26 N=l, NSPEC 

FLUXCC N ) =FLUXCC N ) + DSUMC N ) 

26 CONTINUE 

IFCNEM, LE.NRO+2) GO TO 1 
DO 27 N-l, NSPEC 

I FC C DSUMC N )/FLUXCC N) ) . GT . EPSEM) GO TO 28 

27 CONTINUE 
GO TO 10 

28 CONTINUE 

I CONTINUE 

10 CONTINUE 

DO 31 N=l, NSPEC 

FLUXCC N) =2. DOXXCCN)*FLUXCCN) 

31 CONTINUE 

PRINT 11, NX, N EM, XS, PHI MAX* DEG, QMAX, 

1 C FLUXCC N ), DABS C DL OGIOC FL UXCCN) ) ),N=1, NSPEC) 

II FORMAT C/1X,2I4,F9.4,F7.2,D11.3,5C1X,D10.3, V',F5.2)) 

200 CONTINUE 

PRINT 102 

102 FORMAT C///1X, 'END RINGBD RUN 1 ,///) 

STOP 

END 



RIN0061 
RIN0062 
RIN0063 
RIN0064 
RIN0065 
RIN0066 
RIN0067 
RIN0068 
RIN0069 
RIN0070 
RING07 1 
RIN0072 
RIN0073 
RIN0074 
RIN0075 
RIN0076 
RIN0077 
RIN0078 
RIN0079 
RIN0080 
RIN0081 
RIN0082 
RIN0083 
RIN0084 
R I N 0 0 8 5 
RIN0086 
RIN0087 
RIN0088 
RIN0089 
RIN0090 
RIN0091 
RIN0092 
RIN0093 
RIN0094 
RIN0095 
RIN0096 
RIN0097 
RIN0098 
RIN0099 
RIN0100 
RIN0101 
RIN0102 
RIN0103 
RIN0104 
RIN0105 
RIN0106 
RIN0107 
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SUBROUTINE INIDAT 
IMPLICIT REAL^8CA-H,0-Z,$) 

REALX8 L AMDAO , L AMDA1 

COMMON /GAMA/G,G1 ,G2,G3,G4,G5,G6,G7,G8,G9,G10,G11,G12,G13,G1 



G16,G17,G18,G19,G20 
COMMON /PAR/CO, ENO, EMI ,D,TLIM, ETAL IM , CL IM , EL 0 , QO , TO , 

PBI RD, RBI RD, DM0, DEG, OMEGA, XSVC 51) 

/NPAR/NET A , NC , NT , NEMO , NPHI , NXS , NRO , NSPEC 
/GEOM/APF, PAI,PAI2,SW,CW,BETA,SBETA,CBETA,PSI1,SPSI1, 
CPSI1,PSIF,SPSIF,CPSIF, AK,SK,CK,AO,RF,XF,YF,ZF, 
PHI,SPHI,CPHI,RMIN, RMAX, XS , DIST , 
AMU1,2ETA1,XN,YN,ZN,PSIM,SPSIM,CPSIM,R0 



COMMON 

COMMON 



RIN0108 
RIN0109 
RIN0110 
4, G15, RINO 1 1 1 
RIN0112 
RIN0113 
RIN0114 
RIN0115 
RIN0116 
RINO 1 17 
RIN0118 
RIN0119 
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95 

96 

' 97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 



on 



COMMON /EPSIl/EPSQ, EPSETA , EPST , EPSC , EP5EM 

COMMON /EXTREM/TEXT, ETAEXT , CEXT, REXT, PSI EXT, EMEXT , BEXT , QEXT 
COMMON /SPEC/WAV,XC(5),WC(5),WCR(5),XNAME(5),QC(5),FLUXC(5) 
DATA XC/.091D0, .09 IDO, . 104D0, .135D0, .57 9 DO/ 

DATA WC/1 .00D0,20.0D0,2.0 0D0,21 .0 DO, 4. 0 0 DO/ 

DATA XNAME/' H HF f , f H2 DF f , f HE '/ 

PAI=4.*DATAN( .ID 1) 

AR=8 . 3143D3 
AV = 6 . 022D 26 

C 0MEGA=0 . 5 IS FOR HARD SPHERE COLLISIONS, 

C AN AVERAGE RECOMMENDED VALUE IS ABOUT OMEGA=0.75 
OMEGA=0 . 5D0 
NSPEC=5 
WAV=0. 

DO 51 N = 1 , NSPEC 
WAV=WAV+XC(N)*WC(N) 

51 CONTINUE 

DO 52 N=l, NSPEC 
WCR(N)=DSQRTCWCCN)/WAV) 

52 CONTINUE 
AO = 2 . 5D0 
EMI =4 . 0D0 
RHO0 = 0 . 007 5D0 
T0=2 . 300D3 
D=2 . 50D-1 0 
G=1 . 54D0 

. ENO=RHOO*AV/WAV 

CO=DSQRT(G*AR*TO/WAV) 

PBIRD=0 . 05D0X2.D0 

C RO IS THE RADIUS FOR BEGINNING THE INTEGRATION ALONG THE M=M1 
C CHARACTERI STICC THE AUGMENTED BREAKDOWN SURFACE). 

C NRO IS THE NUMBER OF INTEGRATION INTERVALS ON THIS SEGMENT. 

C FOR NO INTEGRATION ALONG M=M1 CHARACTERISTIC, SET NR0=0. 

RO = 0 . 

NR0=10 
DM0=0 . IDO 

NEM0=20 .D0/DM0+NR0 

TO GET FLUX DUE TO AUGMENTED BREAKDOWN SURFACE SOLELY, ACTIVATE: 
NEM0=NR0 
NPHI =1 0 
NXS=13 
XSI =1 . D-2 
XSF=1 . D1 
XSV(1)=XSI 

IF(NXS.EQ.l) GO TO 111 
DXL=(DLOG(XSF)-DLOG(XSI) )/( NXS-1 . DO ) 

XL I = DLOG( XSI ) 

DO 11 NX=2 , NXS 

XSV(NX)=DEXP(XLI+( NX-1 . DO ) *DXL ) 

II CONTINUE 

III CONTINUE 
EPSEM=l.D-5 
DEG = 18 0 . DO/PAI 
PAI 2=PAI/2 . DO 
GAMMA=G 

G1 =( G-l . DO )/2 . DO 

G2=(G+1 .DO)/ (2. DO* (G-l. DO)) 

G3=G/2 . DO 

G4=(G+1 . DO )/( G-l . DO ) 

G5=DSQRT ( (G+l . DO)/ (G-l . DO) ) 



RIN01200 
RIN01210 
RI NO 1 220 
RI NO 1 230 
RI NO 124 0 
RI NO 1 250 
RINO 1260 
RI NO 127 0 
RI NO 1 28 0 
RI NO 1290 
RI NO-1 30 0 
RIN01310 
RI NO 1 320 
RI NO 1 330 
RI NO 1 340 
RI NO 1350 
RINO 1360 
RI NO 1 37 0 
RI NO 1 38 0 
RIN 01390 
RIN 01400 
RINO 1 4 1 0 
RIN01420 
RIN01430 
RI NO 1 440 
RI NO 1 450 
RINO 1 46 0 
RINO 1 47 0 
RI NO 1 48 0 
RIN01490 
RINO 1500 
RIN0151 0 
RI NO 1 520 
RI NO 1 530 
RI NO 1 54 0 
RINO 1 550 
RI NO 156 0 
RI NO 1 57 0 
RIN 01580 
RIN01590 
RI NO 16 00 
RI NO 16 1 0 
R I N 0 1 6 2 0 
RIN01630 
RIN 01640 
RIN0165Q 
RIN0166 0 
RIN01670 
RIN01680 
RIN01690 
RIN 01700 
RI NO 17 1 0 
RIN 01720 
RIN01730 
RIN01740 
RIN01750 
RIN01760 
RI NO 177 0 
RIN01780 
RIN01790 
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147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 

181 

182 

183 

184 

185 

186 

187 

188 

189 



G6 = 1 .D0/CG-1 .DO) RINQ1J 

G7 =2 . DO/ ( G+l . DO ) RIN01J 

G8=(5.D-1*(G+1 .D0)**2/(G“1 ,DQ))XX(1 . DO/ (G+l . DO)) x RIUOU 

1 ( (G+l .D0)/(G-1 . DO))XX( (G-l .DO)/ (G+l .DO)) RINOl^ 

G9=(G+3.D0)/(2.D0*(G“1 .DO)) RINOli 

G10=(7.D0-3.D0*G)/(2.D0*(G“1.D0)) RIN01J 

G11=DSQRT(G/PAI)/(2.D0X(G+1 .DO)) RIN01J 

G12=DSQRT(G/2.D0) RINOl^ 

G13=l . DO/DSQRT ( 2 . D0XGXPAIXX3) RINOU 

LAMDA0=1 ,D0/(DSQRT(2 .D0)XPAIXDXX2XEN0) RINOU 

LAMDA1=LAMDA0*(1 .D0+G1*EM1*X2)XX(G6-0MEGA+0.5D0) RINOU 

RBIRD-Gll/C D*X2*EN0XPBI RD) RINOU 

ZETA1=G5*DATAN(DSQRT(EM1X*2-1 .DO )/G5) RINOU 

AMU1 = DARSI N( 1 .DO/ EMI) RINOU 

PSI 1 =PAI2+AMU1 RINOU 

SPSI1=DSIN(PSI1) * RINOU 

CPS 1 1 =DCOS( PSI 1 ) RINOU 

PSIF=PAI2+AMU1+ZETA1-G5XPAI2 RINOU 

SPSIF = DSIN(PSIF) RINOU 

CPSIF=DCOS(PSIF) RINOU 

CALL BREAKR(EM1,RMIN) RIN02I 

RSMIN=RMIN*DSQRT( (2. DO/ (G-l . DO )+EMl*X2)/( EM1XX2-1 -DO)) RIN02( 

C RIN02I 

PRINT 201,N$PEC,XNAME RIN02I 

201 F0RMATC/1X, ’SPECIES DATA NSPEC=’,I3/ RIN02I 

1 IX, ’SPECIES NAMES * , 1 1 ( 2X , A6 , 2X ) ) RIN02I 

PRINT 202, XC RIN02I 

202 FORMAT ( IX, ’MOLE FRACTION XC= SI 1 ( F8 . 4 , 2X) ) RIN02I 

PRINT 203, WC RIN02I 

203 FORMAT ( IX, ’MOL. WEIGHT WC= ’ , 1 1 ( F8 . 4 , 2X ) ) RIN02I 

PRINT 21,AR,AV,WAV,G,RHOO,TO,ENO,CO,D RIN02! 

21 F0RMATC/1X, ’THERMODYNAMIC DATA’/ RIN02: 

1 IX, ’ AR, AV , WAV, GAMMA = ’,2X,2D14.5,2F9 . 3/ RIN02! 

2 IX, ’RHOO,TO, ENO,CO, D= ’ ,D12.4, F8 . 0 , D1 3 . 5 , 2D1 2 . 4 ) RIN02! 

PRINT 22,EM1,PSI1*DEG,PSIF*DEG,PBIRD, RIN02! 

1 AO,RMIN,RSMIN,RO, RIN02! 

2 LAMDAO , LAMDA1 , RMIN/LAMDA1 , RSMIN/L AMDA1 RIN02! 

22 F0RMATC/1X, ’FLOW AND GEOMETRY DATA’/ RIN02! 

1 IX, ’EM1,PSI1,PSIF,PBIRD=’,4F9.3/ RIN02: 

2 IX, ’A0,RMIN,RSMIN,R0=’,F9.3,3D14.5/ RIN02! 

3 IX, ’ LAMDAO, LAMDA1 , RMIN/L AMDA1 , RSMI N/ L AMDA1 = ’ , 4D13 . 4) RIN02; 

PRINT 23, DMO,NPHI ,NRO, EPS EM RIN02; 

23 F0RMATC/1X, ’INTEGRATION DATA’/ RIN02; 

1 IX, ’ DMO,NPHI,NRO, EPSEM=* ,F9 .3,215, 4X, D12.3) RIN02J 

RETURN RIN02; 

END RIN022 

SUBROUTINE FLUX( EM, R ) RIN02; 

IMPLICIT REAL*8(A-H,0-Z,$) RIN02; 



COMMON /GAMA/G, G1 ,G2,G3,G4,G5,G6,G7,G8,G9,G10,G11,G12,G13,G14,G15,RIN02< 

RIN022 
RIN02: 
RIN02: 
RIN02: 
RIN02: 
RIN02: 
RIN02: 
RIN02: 
RIN02: 
RIN02: 



1 G16,G17,G18,G19,G20 

COMMON /PAR/CO, ENO, EMi, D,TLIM, ETAL IM, CLIM, ELO , QO, TO , 

1 PBIRD, RBIRD, DMO , DEG, OMEGA , XSV( 51 ) 

COMMON /NPAR/NETA , NC, NT, NEMO , NPHI , NXS , NRO , NSPEC 

COMMON /GEOM/APF,PAI ,PAI2,SW,CW, B ETA , SB ETA , CBETA , PSI 1 , SPSI 1 , 

1 CPSI1,PSIF,SPSIF,CPSIF,AK,SK,CK,A0,RF,XF,YF,ZF, 

2 PHI , SPHI , CPHI , RMIN, RMAX, XS, DIST, 

3 AMU1,ZETA1,XN,YN,ZN,PSIM,SPSIM,CPSIM, RO 
COMMON /EXTREM/TEXT, ETAEXT,CEXT, REXT, PS I EXT , EMEXT , BEXT , QEXT 
COMMON /SPEC/ WAV , XC(5),WC(5),WCR(5), XNAMEC 5 ) , QCC 5 ) , FLUXCC 5 ) 
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190 

191 

192 

193 

194 

195 

196 

197 

193 

199 

200 

201 

202 

203 

204 

205 

206 

207 

203 

209 

210 

211 

212 

213 

214 

215 

216 

217 

218 

219 

220 

221 

222 

223 

224 

225 

226 

227 

223 

229 

230 

231 

232 

233 

234 

235 

236 



SPHI = DSI N ( PHI ) 

CPHI = DCOS ( PHI ) 

XMEAN=R*CPSIM 
YMEAN=(AO+R*SPSIM)*CPHI 
ZMEAN=(AO+R*SPSIM)*SPHI 
T BET A =ZMEAN/ (YMEAN-AO ) 

BETA=PAI2-DATAN( 1 . DO/ T BET A) 

SBETA=DSIN(BETA) 

CBETA=DCOS(BETA) 

DIST=DSQRT( ( XS-XMEAN)**2+( YMEAN-AO ) **2+ZMEAN**2 ) 

CW=-(XS-XMEAN)/DIST 

SW=DSQRT( 1 . D0-CWXX2 ) 

GM = C 1 . D0+G1*EMXX2)X*(-G2) 

AMU=DARSINC1 .DO/EM) 

TETA=PSIM-AMU 
STETA=DSIN( TETA ) 

CT ETA = DCOS ( T ETA ) 

CKAPA=(CTETA)X(-CW)+(STETA*CPHI)*(-SWXCBETA)+ 

1 (STETAXSPHI)X(-SWXSBETA) 

SKAPA=DSQRT(1 . D0-CKAPAXX2 ) 

QEXT = 0 . 

DO 1 N=1 , NSPEC 
EMT=EM*CKAPA*G12XWCR(N) 

I F( DABSC EMT ) . GT . 13 . DO ) EMT = EMT* ( 1 3 . DO/DABS ( EMT ) ) 

P0N=G3*EM*X2*WCR(N)XX2 

P0NT=P0W*SKAPAXX2 

IFC POW.GT . 1 .0D2)P0W = 1 .0D2 

IFCPOWT.GT.l .0D2) POWT=1.0D2 

EXP1 =DEXP ( -POW ) 

EXP2=DEXP ( -PONT ) 

ERFC1=DERFC( -EMT ) 

IFCERFCl.LT. l.D-43) ERFCl=l.D-43 
C IFCXS.LT. 5. DO) GO TO 234 

C IFCEXPl .GT. 1 . D-20 . AND.ERFC1 .GT.l .D-20) GO TO 234 

C PRINT 235 , N , PON, POWT , EMT , EXP1 , EXP2, ERFC1 

C235 F0RMAT(I4,3D11 .3) 

C234 CONTINUE 

EVER1=EN0X(C0/WCR(N))XG13XGMX(1 .D0+EMTXX2)XEXP1 
EVER2=( ENOXCCO/WCRCN) )xo. 5D0/PAI)X(EMXWCR(N) )XCKAPAXGMX 
1 (1 .5D0+EMTXX2)XEXP2XERFC1 

QCC N) = EVER1 + EVER2 
I F( QEXT . GE . QCC N ) ) GO TO 1 
QEXT=QC(N) 

1 CONTINUE 

RETURN 
END 



RIN02390 
RIN 02400 
RIN 02410 
RIN02420 
RIN02430 
RIN02440 
RIN 02450 
RIN02460 
RIN 02470 
RIN 02480 
RIN02490 
RIN02500 
RIN0251 0 
RIN02520 
RIN02530 
RIN 02540 
RIN02550 
RIN02560 
RIN 0257 0 
RIN02530 
RIN02590 
RIN02600 
RIN0261Q 
RIN02620 
RIN02630 
RIN02640 
RIN02650 
RIN02660 
RIN02670 
RIN02680 
RI NO 26 9 0 
RIN027 00 
RI N 027 1 0 
RIN027 20 
RIN02730 
RIN 027 40 
RIN027 50 
RIN027 6 0 
RIN0277 0 
RIN02730 
RIN02790 
RIN02300 
RIN02810 
RIN02320 
RIN02830 
RIN02840 



SUBROUTINE BREAKRC EM, R ) 

IMPLICIT REALX3(A-H,0-Z,$) 

COMMON /GEOM/APF , PAI , PAI 2 , SW , CW , BETA , SBETA , CB ETA , PSI 1 , SPSI 1 , 

1 CPSI1,PSIF,SPSIF,CPSIF,AK,SK,CK, AO , RF, XF, YF, ZF , 

2 PHI, SPHi, CPHI, RMIN, RMAX , XS , DI ST , 

3 AMU1,ZETA1,XN,YN,ZN,PSIM,SPSIM,CPSIM,R0 
COMMON /PAR/CO, ENO, EMI, D,TLIM, ETAL IM, CLIM, ELO,QO,TO, 

1 PBIRD, RBI RD, DMO, DEG, OMEGA, XSVC 51) 

COMMON /GAMA/G,G1,G2,G3,G4,G5,G6,G7,G3,G9,G10,G11,G12,G13,G14 
1 G16 , G17 ,G18,G19,G20 

R = RBIRDXDSQRT( EMXX2-1 .DO)X( 1 . DO + Glx EMX)(2 ) ( G6 -OMEGA + O . 5D0 ) 

ZETA=G5*DATAN(DSQRTCEM*X2-1 .D0)/G5) 

PSI=PAI2+AMU1+ZETA1-ZETA 



RIN02850 
RIN028 6 0 
RIN028 7 0 
RIN 02880 
RIN0289 0 
RIN 02900 
RIN0291 0 
RIN02920 
, G15 , RIN02930 
RIN02940 
RIN 02950 
RIN0296 0 
RIN0297 0 
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238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 



XF=R*DCOS(PSI) 

YF=R*DSIN(PSI)+A0 
ZF = 0 . 

1 CONTINUE 
RETURN 

ENTRY BREAKM(EM,R) 

R=RBIRD*DSQRT(EMXX2-1 .D0)x( 1 .D0+Gl*EM*X2)*x(G6~OMEGA+0.5D0) 
ZETA=G5 xDATAN(DSQRTCEM*X2-1 .D0)/G5) 

PSIM=PAI2+AMU1+ZETA1-ZETA 

SPSIM=DSINCPSIM) 

CPSIM=DCOS(PSIM) 

RETURN 

END 



RI NO 298 0 
RI NO 29 9 0 I 
RIN03000 : 
RIN 03010 
RIN03020 
RIN 03030 
RIN03040 
R I N 0 3 0 5 0 
RIN03060 
RIN0307 0 
RI NO 3 08 0 
RING 30 9 0 
RIN031 00 
RIN 03110 
RIN03120 



GENTRY 



RIN031 30 



- FLUX INTEGRATION FROM BREAKDOWN SURFACE 
SPECIES DATA NSPEC” 5 



SPECIES NAMES 
MOLE FRACTION 
MOL. WEIGHT 


r x 
nn 

ii h i 


H 

0.0910 

1.0000 


HF H2 

0.0910 0.1040 

20.0000 2.0000 


OF 

0.1350 

21.0000 


HE 

0.5790 

4.0000 


THERMODYNAMIC DATA 
AR.AV.WAV, GAMMA- 


0.831430 04 0 602200 27 


7.270 


1.540 



RHDO . TO , ENO , CO , D 3 0.7500D-02 2300, 0.621250 24 0.20130 04 0.25000-09 

FLOW AND GEOMETRY DATA 

EMI , PS 1 1 ,PS1F,PB1RD= 4.000 104.478 41.044 0.100 

A0.RM1N, RSM1N, RO 2 2.500 0.303740-02 0.348120-02 0.000000 00 

*.AMDA0*LAMDA1,RMIN/LAMDA1,RSMIN/LAMDA1- 0.57970-05 0.12810-03 0.23720 02 0.27180 02 

INTEGRATION DATA 

DMO* NPHl * NRO , EPS EM 3 0.100 10 10 0.1000-04 



NX 


NEM 


XS 


PHI MAX 


QMAX 




H 


/ LOG 


HF 


/ 


LOG 


H2 


/ LOG 


OF 


/ 


LOG 


HE 


/ LOG 


1 


47 


0.0100 


11.19 


0 .1200 


25 


0.8240 


22/21.92 


0.2690 


07/ 


6.43 


0 .777 0 


21/20.89 


0.6910 


06/ 


5.84 


0.6530 


20/19.81 


2 


50 


0.0178 


12.15 


0.1210 


25 


0 . 4740 


22/21.68 


0.1620 


07/ 


6.21 


0.4240 


21/20.63 


0.4180 


06/ 


5.62 


0.3490 


20/19.54 


3 


54 


0 0316 


13.50 


0.1220 


25 


0.2370 


22/21.38 


0.7760 


06/ 


5.89 


0.1960 


21/20.29 


0.2020 


06/ 


5.30 


0.1540 


20/19.19 


4 


59 


0 . 0562 


15.27 


0 . 1230 


25 


0.1060 


22/21.02 


0.2950 


06/ 


5.47 


0.7800 


20/19.89 


0.7700 


05/ 


4.89 


0. 5680 


19/18.75 


5 


64 


0 .1000 


17 . 16 


0. 1240 


25 


0.4210 


21/20.62 


0.8920 


0 5/ 


4.95 


0.2680 


20/19.43 


0.2330 


0 5/ 


4.37 


0. 1750 


19/18.24 


6 


70 


0.1778 


19.55 


0 . 1240 


25 


0.1490 


21/20.17 


0.2190 


0 5/ 


4.34 


0.7910 


19/18.90 


0. 5750 


0 4/ 


3.76 


0.4510 


18/17 .65 


7 


77 


0.3162 


22.52 


0.1240 


25 


0.4710 


20/19 67 


0 . 4640 


04/ 


3.67 


0.2030 


19/18.31 


0 . 1220 


04/ 


3.09 


0 .1000 


18/17.00 


8 


84 


0.5623 


25.65 


0.1240 


25 


0.1320 


20/19.12 


0.9010 


03/ 


2.95 


0.4610 


18/17 .66 


0.2370 


03/ 


2.37 


0.2010 


17/16.30 


9 


92 


1.0000 


29 . 39 


0. 124D 


25 


0.3350 


19/18.52 


0.1670 


0 3/ 


2.22 


0.9590 


17/16.98 


0.4400 


02/ 


1.64 


0.3810 


16/15.58 


10 


100 


1.7783 


33.26 


0.1240 


25 


0.7710 


18/17.89 


0 .3050 


02/ 


1.48 


0 1870 


17/16.27 


0.8010 


01/ 


0.90 


0.6990 


15/14.84 


11 


108 


3.1623 


37.22 


0.1240 


25 


0.1650 


18/17.22 


0.5480 


01/ 


0.74 


0.3530 


16/15.55 


0.1440 


01/ 


0.16 


0.1260 


15/14.10 


12 


116 


5.6234 


41 . 19 


0 . 124D 


25 


0.3320 


17/16 52 


0.9810 


00/ 


0 .01 


0.6480 


15/14.81 


0.2580 


00/ 


0.59 


0 .227D 


14/13.36 


13 


124 


10 . 0000 


45.14 


0 .1240 


25 


0.6410 


16/15.81 


0 1750 


00/ 


0.76 


0.1170 


15/14.07 


0 .4610-01/ 


1.34 


0 .405D 


13/12.61 



END R1NGBD RUN 
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Figure 1. Thermal Backseat tering from Laser Exhaust Plume 
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Figure 2. Prandtl-Meyer Centered Rarefaction 



Fan and 



Breakdown Surface 



( Schematic) 
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Figure 3. Flux Integration Scheme 
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STREAMLINE 2 



M 



Figure 4 



I 



R(M 1 )=54.4*A 1 
-BREAKDOWN SURFACE R(M ] )=47.4*A 1 

'STREAMLINE 1 R(M ] )=47.4*A 1 




Prand tl-Meyer Flow Field Near the Corner, Including Actual 

Streamlines and Breakdown Surface 
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FLUX OF SPECIES H 
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1023 , 



M-| =4 y-1.54 

Tq = 2300 (K) W a = 7.27 B = .05 




Figure 5. Flux of Species H at Various Stagnation Densities 
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FLUX OF SPECIES HF 
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Figure 6. Flux of Species HF at Various Stagnation Densities 
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Figure 7. Flux of Species H 2 at Various Stagnation Densities 
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Figure 8. Flux of Species DF at Various Stagnation Densities 
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Figure 9, Flux of Species He at Various Stagnation Densities 
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Figure 10. Flux of Species HF at Various Values of Breakdown Parameter 
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Figure 11. Flux of All Species at Typical Operating Conditions 
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Figure 12. Schematic Display of Complete Breakdown Surface in 
Exhaus t Plume 



Ringjet 
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